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ABSTRACT 

Solar tri-diurnal variations of cosmic rays have been 
analyzed in a wide range of rigidity, using data from neutron 
monitors, and the surface and underground muon telescopes for 
the period 1978-1983. The rigidity spectrum of the aniso- 
tropy in space is assumed to be of power-exponential type as 

(P/Y p o)^exp(y-P/P 0 ) . By means of the best-fit method between 
the observed and the expected variations, it is obtained that 
the spectrum has a peak at P (=yP 0 )=9O GV, where y=3.0 and 
Po=30 GV. The phase in space of the tri-diurnal variation 
is also obtained as 7.0 hr (15 hr and 23 hr LT) , which is 
quite different from that of M. hr, arising from the axis- 
symmetric distribution of cosmic rays with respect to the IMF. 

1. Introduction 

Studies of higher harmonics in the daily variation of cosmic rays 
provide valuable infromations as to the nature of the cosmic ray modula- 
tion in the heliosphere. So far a great many investigation has been 
performed of the solar tri-diurnal variation, and some of its character- 
istic natures have been revealed; its extra-terrestrial origin and also 
its rising spectrum with rigidity (e.g., Bieber et al., 1983; Mori et al., 
1984; hereafter referred to as Paper I). 

These results, however, have not yet been firmly founded because of 
the statistical uncertainty of the data used which is largely due to 
limited rigidity range. In particular, there is a wide divergence of 
opinion as to the direction of the anisotropy; some analyses reported 
*\,1 hr direction, while some obtained %7 hr direction. This discrepancy 
is crucial for determining the modulation mechanism; 1 hr direction is in 
favor of the * loss-cone* model (Fujii, 1971) or the model of Bieber et al. 
(1983) , whereas 7 hr direction requires other kind of model such as the 
one presented by Munakata and Nagashima (1985; in this issue). 

In the present analysis, we try to obtain a definite form of the 
anisotropy responsible for the averaged tri-diurnal variations over 
6 years (1978-1985) , using data in a wide range of rigidity from neutron 
monitors to underground muon telescopes, which respond to the effective 
rigidity of 10 to 600 GV. 

2. Analysis of Data 

The present experimental data include 2 to 6-year averages of neutron 
monitor data in a worldwide network (36 stations from the polar to the 
equatorial regions) and 6-year averages of the surface (NAGOYA) and under- 
ground (MISAT0 and SAKASHITA) muon telescopes for the period 1978-1983. 
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Pig. 1 illustrates the summation harmonic. dial, showing clearly the 
persistency of year-to-year changes of the observed tri-diurnal variations 
of three muon telescopes at NAGOYA (surface) for 1971-1983, at MISATO 

(34 m.w.e. in depth underground) 
for 1974-1983 and at SAKASHITA 

SOLAR TRI-DIURNAL variation (80 m.w.e. in depth underground) 

(Vertical component) ^ for 1978-1983. In the present 

analysis, 6-year averages over 
o o.o s% — jr- 2 1978-1983 were used; of 13 direc- 

^ y l tional components at NAGOYA, 9 com- 

0 h ponents at MISATO and 10 components 

\ y at SAKASHITA, and these are listed 

ItT in Paper I . 

/ 74 

Hisato s' 1 * i In orc ^ er t 0 derive the 3rd 

— S j _ \ order anisotropy, the tri-diurnal 

/ variation produced by the aniso- 

3 01 I tropy with a certain form of the 

Sakashita / rigidity spectrum G(P), will be 

^ 8l compared with the observed results. 

Nagoya In the anaysis, G(P) is assumed to 

be of power-exponential type as 

Fig, 1 Observed solar tri- 

diurnal variations of V- G(P)=(P/yPo) exp(y-P/Po) (1) 

component at NAGOYA, MISATO 

and SAKASHITA in summation Evidently, G(P) becomes maximum at 

dial. P=yPo for y>0, and is normalized 

to unity. 


Fig, 1 Observed solar tri- 
diurnal variations of V- 
component at NAGOYA, MISATO 
and SAKASHITA in summation 
dial. 


The expected 3rd harmonic components (A. ., B, . ) for j-th components 

1 1 i D 

telescopes of i-th station including neutron monitor station, are related 
to the space harmonic components (X, Y) by the following equations with 
coupling coefficients c, . and s. . as 


A. . - (c. .X + s. .Y) # 
1 3 13 ID 


(— s * »X + c. >Y) 
ID ID 


The coefficients and were derived from the differential coupling 

coefficients dc, . and ds . . given by Fujimoto et al. (1984). The best-fit 
ID ID 

parameters satisfying the following equation are determined. 

E.w. .{ (A. . - a. ,)^ + (B. . - b. ,)^ }- minimum = 6^ (3) 

ID ID iD ID ID ID 

In Eq. (3) , the weight w. , is given for the observed a, . and b. . so as to 

ID ID ID 

balance the difference in the number of the data used in the rigidity 
interval; 6 divisions in all the ranges of 10 to 600 GV to have equal 
weight. Also each vector in all the muon components was subtracted from 
the reference component to eliminate the atmospheric temperatute effect; 
in this analysis, V-component for NAGOYA and MISATO and SE-component for 
SAKASHITA were taken as the reference. 


3. Results and Discussions 
■ — 2 

The contour map of equal 6 *s obtained from the residual sum of the 
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observed vectors and the expected vectors computed by Eq. (2), is dhown 


against y and P 0 in Fig. 2, 



As is seen in the figure, minimal <$ region 
locates in relatively narrow 
regions of y and P 0 -value. It is 
noted that in these regions the 
product of y and Pq keeps almost 
invariant as yPo=90 GV. The best- 
fit parameters were obtained as 
that y=3. 0 and P 0 = 3O GV. 

The direction (phase) of the 
anisotropy was obatined as 7.0 hr. 
It is also noted that this direc- 
tion shows almost constant value 
for the above minimal region. 


The fitness of the reproduced 
variations computed by Eq. (2) in 
the best-fit case to the observed 
ones is shown in Fig. 3; (a) for 

NAGOYA , (b) for MI SATO and (c) for 

SAKASHITA, respectively. The 
observed vectors are shown with 
the solid circles (•) and the 

reproduced ones with the open circles (O). Also the reference component 
mentioned above is shown by the symbol (©) and the origin of the diagram 
in the best-fit case is given by the cross (X). A fairly cjood agreement 
can be seen between two patterns drawn by connecting each observed point 
(solid line) and each reproduced point (dotted line) for each station. 


Fig, 2 Contour map of equal 6* 
plotted against y and P , 
showing minimal <$ 2 (=8) for 
y=3 . 0 and P 0 s3O GV. 


The present anaysis can obtain a definite parameter set as; y=3.0 
and Po=30 GV. As shown in Paper I , where only muon data at the above 
three stations (NAGOYA , MISATO and SAKASHITA) were used, parameter set 
(Y/Po) rather spereads as; y^2^12 and Po^SO^S GV. The present form can be 
derived by adding neutron monitor data to Paper I . Note that the observed 
vectors of neutron monitors from 36 stations are so divergent even if in 
the best-fit case, which may be largely due to inaccurate barometric cor- 
rection in the data. 


On the other hand, the direction (phase) of the anisotropy in space 
was obtained as ^7 hr (VL5 hr and <v23 hr LT) . So far there has been a 
wide divergence of opinion as to the direction; 'Vl hr direction was 
reported by the analyses of neutron monitor data alone, while V7 hr 
direction was reported by using muon data. The former is in favor of the 
’loss-cone 1 model (Fujii, 1971) or the model of Bieber et al. (1983), 
arising from the axis-symmetric distribution of cosmic rays with respect 
to the IMF. On the other hand, the latter requires other kind of model 
such as the one preseneted by Munakata and Nagashima (1985, in this issue). 
The present result provides a firm basis for the latter model. 

A further analysis would be necessary to derive the detailed form and 
its year-to-year variation. 
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(a) for NAGOYA 


0 H 

| % SAKASHITA 

t 0.02 



0 H 



Fig. 3 Fitness between the 
observed ( 0 ) and the re« 
produced (O) variations; 

(a) for NAGOYA, (b) for MI SATO 
and (c) fro SAKASHITA for the 
period 1978-1983* 


(c) for SAKASHITA 
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